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Abstract 

The unbound nucleus ^^Li is evaluated by studying three-neutron one-proton 
excitations within the multistep shell model in the complex energy plane. It 
is found that the ground state of this system consists of an antibound 2~ 
state. A number of narrow states at low energy are found which ensue from 
the coupling of resonances in ^^Li to continuum states close to threshold. 
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1. Introduction 

The study of halo nuclei is one of the main subjects of research in nu- 
clear physics at present. Many theoretical predictions on halo, superhalo 
and antihalo nuclei have been advanced in recent years [H, Q, IH. Most of 
these calculations correspond to nuclei very far from the stability line. They 
are mainly thought as a guide for experiments to be performed in coming 
facilities. The general feature found in these calculations is that a necessary 
condition for a nucleus to develop a halo is that the outmost nucleons move 
in shells which extend far in space. That is, only a weak barrier keep the sys- 
tem within the nuclear volume. These shells may be resonances, antibound 
states (also called virtual states), or even low-spin bound states which lie 
very close to the continuum threshold. These conditions are fulfilled by the 
nucleus ^^Li and also heavier Li isotopes. There are a number of experiments 
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which have been performed in these very unstable isotopes in order to get 
information about the structure of halos [4] • In particular we will concentrate 
our attention to Refs. jsi 0, 0] where the spectrum of ^^Li was measured. 
Our aim is to analyze these experimental data by using a suitable formalism 
to treat unstable nuclei. This formalism is an extension of the shell model 
to the complex energy plane and is therefore called cornplex shell model jsj, 
although the name Gamow shell model is also used [9[. In addition, the 
correlations induced by the pairing force acting upon particles moving in de- 
caying single-particle states will be taken into account by using the multistep 
shell model (MSM) fu^. 

The formalism is presented in Section |21 Applications are in Section [3] 
and a summary and conclusions are in Section HI 



2. The formalism 

The study of unstable nuclei is a very difficult undertaking since, in prin- 
ciple, time dependent formalisms should be used to describe the motion of 
a decaying nucleus. However, the system may be considered stationary if it 
lives a long time. In this case the time dependence can be circumvented. In 
fact, often unstable nuclei live a very long time and therefore they may be 
considered bound as, e.g., in alpha decaying states of many heavy isotopes, 
like ^°^Bi or ^^''Ta(9~), with T1/2 > lO^^y. On the other hand, experimental 
facilities allow one nowadays to measure systems living a very short time. To 
describe these short time processes one has to consider the decaying character 
of the system. 

Of the various theories that have been conceived to analyze unbound 
systems, we will apply an extension of the shell model to the complex energy 
plane |8j]. The basic assumption of this theory is that resonances can be 
described in terms of states lying in the complex energy plane. The real parts 
of the corresponding energies are the positions of the resonances while the 
imaginary parts are minus twice the corresponding widths, as it was proposed 



by Gamow at the beginning of quantum mechanics These complex states 
correspond to solutions of the Schrodinger equation with outgoing boundary 
conditions. We will not present here the formalism in detail, since this was 



done many times before, e.g., in Refs. [12|, ll3|. Rather, we will give the main 



points necessary for the presentation of the applications. 
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2.1. The Berggren representation 

In this Subsection we will very briefly describe the representation to be 
used here. 

The eigenstates of a central potential obtained as outgoing solutions of 



the Schrodinger equation can be used to express the Dirac 5-function as [14 
A(r-/)^5:u,„(r)u,„(/)+/ .^.(r.^).(r'.E). (1) 

where the sum runs over all the bound and antibound states plus the complex 
states (resonances) which lie between the real energy axis and the integration 
contour L~^. The wave function of a state n in these discrete set is w„(r) 
and u{r, E) is the scattering function at energy E. The antibound states 
are virtual states with negative scattering length. They are fundamental to 
describe nuclei in the Li region jl5| . 

The resonances and the antibound states are poles of the single-particle 
Green function and, therefore, we will call them "poles" in order to make a 
distinction with the scattering states. This is important since it is through 
the poles that we will recognize physically meaningful states, i. e. states that 
live a time long enough. Although meaningful states are usually immersed 
among continuum states, they can be recognized because their wave func- 
tions contain important contributions from the poles, as will be seen in the 
Applications. 

Discretizing the integral of Eq. ([1]) one obtains the set of orthonormal 



vectors \ipj) forming the Berggren representation [16j. Since this discretiza- 
tion provides an approximate value of the integral, the Berggren vectors fulflU 
the relation I ^ \ipj){ipj\, where all states, that is bound, antibound, res- 
onances and discretized scattering states, are included. The corresponding 
single-particle wave functions are 

(r>,) = i?„,,,,,(r)(xi/2n(0),,„. (2) 
where x is the spin wave function and 

Rn,kki.'r) = (i)n,Uki.'r)/r (3) 

is the radial wave function fulflUing the Berggren metric, according to which 
the scalar product between two functions consists of one function times the 
other (for details see Ref. [l^), i.e., 

dr(i)n,Uj, {r)(pn(kn (^) = ■ (4) 
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We will apply the Berggen representation to analyze the spectrum of ^^Li 
by using the Multistep Shell Model Method [lo| (MSM). In order to make 
the presentation clear we will give a short description of the main points to 
be used in this paper. 

The Berggren representation has been used before within the framework 
of the CXSM (or the Gamow Shell Model, which is the same [l3[) to study Li 
isotopes (17']. The core in these calculations was assumed to be '^He and the 
Shell Model single-particle states were the shells Op only. As pointed out in 
that reference, the study of many particles moving in states lying in the com- 
plex energy plane can be a challenging task. The main problem is that due 
to the presence of the scattering states the dimension of the Berggren basis 
soon becomes very large, as well as non-Hermitian and complex. Therefore 
in 171] one could describe well light Li isotopes, up to ^Li. But in heavier 



isotopes the state lsi/2 (which is an antibound state) is of a fundamental 
importance and its inclusion would make the application of the CXSM a 
prohibitive undertaking The importance of the antibound state in this 



nuclear region was known since a rather long time [15|. An excellent expla- 



nation of the halo nucleus ^^Li was given in Ref. 18| by using the Continuum 
Shell Model, including in the representation only the neutron waves si/2, pi/2 
and 0^3/2, which in terms of the CXSM implies to include the antiboud state 
and the resonances Opi/2 and 0^3/2, but not the state Op^/o. The correspond- 



ing CXSM calculation was indeed performed in Refs. [19|, l20[ and, perhaps 
not surprising, one could thus explain well the structure of ^^Li as well as the 
corresponding halo. 

One way of avoiding too large dimensions is by including in the basis only 
physically meaningful states. These are states which govern the calculated 
quantities. For instance, in the evaluation of ^Li the single-particle resonance 
Op\/2 is very broad and its inclusion does not affect the results appreciably 



17|. In this context, one important point that has to be emphasized is 
the evolution of the single-particle states in these isotopes as the continuum 
threshold is reached. In light Li isotopes the unstable character of the un- 
bound states is reflected by the CXSM in that the neutron resonance Opi/2 
is very broad and plays practically no role in the evaluation of the spectrum. 
Instead, starting in ^Li this state becomes extremely important and in ^^Li 



contributes by about 50% to the formation of the halo 20|, |2l| . On the other 



hand, the state 0p3/2 is fundamental in the structure of light Li isotopes, but 
its importance is diminished to a point where it can be neglected starting in 
^°Li. Shell Model calculations performed within this approximation repro- 
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duced very well the experimental data 18|, |20|, |22| . Even three-body (Fadeev) 



approaches 2l|, |23| have shown that the structure of heavy Li isotopes are 
not appreciably affected by the state 0p3/2. 

An important conclusion of all the calculations mentioned above is that 
the low lying states in ^^Li are essentially two-neutron excitations. Moreover, 
the corresponding single-neutron states, forming the CXSM representation, 
are in a similar way extracted from the odd-neutron excitations in ^°Li [2^, 



21| . It is to be noticed that in these references the proton degrees of freedom 
were also neglected. The reason for this is that the pairing interaction acting 
upon the even number of neutrons present in odd Li isotopes determines the 
spectrum, leaving the protons as mere spectators. However, this cannot be 
said in our case of ^^Li. We will therefore include in the calculations the 
proton shell 0^3/2. 

In our three-neutron one-proton case a way of distinguishing physically 
meaningful states is by considering first the two-neutron excitations. From 
the calculated states one can single-out the resonances which are very broad 
or otherwise unphysical (for instance energy eigenvalues with real as well as 
imaginary negative parts). 

In order to introduce the physically meaningful two-neutron states in 
the formalism we will apply the MSM. In this method one solves the Shell 
Model equations in several steps. One first chooses a single-particle repre- 
sentation. Then one evaluates the two-body equations. Next one solves the 
three-body equations within a basis consisting of the tensorial product of the 
one- and two-body basis previously evaluated. For the four-particle case one 
can choose as a basis the one- times three-particle basis states already evalu- 
ated or the two- times two-particle basis, and so on. In our case we will first 
solve the two- neutron states as done in Refs. [2o|. We will thus sort out the 
physical meaningful states to be used in the three-neutron case. Besides the 
advantage of reducing the basis dimensions, the inclusion of these meaningful 
two-particle states implies that relevant continuum (scattering) states, that is 
those continuum states that determine the physically meaningful two-particle 
resonances, are also included. 

In the next MSM step we will evaluate the three-neutron states in a 
basis consisting of the one- times two-particle states. Finally, with the three- 
neutron states thus determined we will form the one-proton three-neutron 
MSM basis to evaluate the spectrum of ^^Li. 

The formalism corresponding to these calculations starts by choosing the 
single-particle (Berggren) states. Using the Berggren representation thus 
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chosen one gets the two-particle shell-model equations in the complex energy 
plane (CXSM) [l2|, i.e., 



{W{a2) - - ej)X{ij; as) = "^{kl; a2\V\ij; a2)X{kl; aa), (5) 

k<l 

where V is the residual interaction. The tilde in the interaction matrix 
element denotes mirror states so that in the corresponding radial integral 
there is not any complex conjugate, as required by the Berggren metric. The 
two-particle states are labeled by 02 and Latin letters label single-particle 
states. W{a2) is the correlated two-particle energy and is single-particle 
energy. The two-particle wave function is given by 

|«2) =P+(«2)|0), (6) 

where the two-particle creation operator is given by, 

P^a2) = J2x{^J;a2)^^, (7) 
and is the angular momentum of the two-particle state. 



We will use a separable interaction, as in Ref. [20[, which describes well 
the two-neutron states in ^^Li. The energies are thus obtained by solving the 
correspondin g di spersion relation. The two-particle wave function amplitudes 
are given by 20 

X{^J■a2) = N^,^^^^ (8) 

where f{ij, ^2) is the single particle matrix element of the field defining the 
separable interaction and is the normalization constant determined by 
the condition ^j<j -^(^j; ^2)^ = 1- 

With the two-neutron states thus evaluated we proceed to the calculation 
of the three-neutron states by using the MSM in the complex energy plane 
(CXMSM), as briefly described below. 

2.2. The Multistep Shell Model Method 

As its name indicates, the Multistep Shell Model Method (MSM) solves 
the shell model equations in several steps. In the first step the single-particle 
representation is chosen. In the second step the energies and wave functions 
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of the two-particle system are evaluated by using a given two-particle interac- 
tion. The three-particle states are evaluated in terms of a basis consisting of 
the tensorial product of the one- and two-particle states previously obtained. 
In this step the interaction does not appear explicitly in the formalism. In- 
stead, it is the wave functions and energies of the components of the MSM 
basis that replace the interaction. The MSM basis is overcomplete and non- 
orthogonal. To correct this one needs to evaluate the overlap matrix among 



the basis states. A general description of the formalism is in Ref. [10||. The 



particular system that is of our interest here, i.e., the three-particle case, can 



be found in Ref. [2J] , where the MSM was applied to study the three- neutron 
hole states in the nucleus ^°^Pb. 

Using the Berggren single-particle representation described above, we will 
evaluate the complex energies and wave functions of ^^Li using the MSM basis 
states consisting of the Berggren one-particular states, which are states in 
^°Li, times the two-neutron excitations that determine the low lying spectrum 
of ^^Li. Below we refer to this formalism as CXMSM. 



The three-particle energies 1^(03) are given by 24 

= E I - - £fc)A(m2, J/32; A;) I (a3|(c;P+(/32)).3|0), (9) 

where Latin letters label single-particle states and e are the corresponding 
single-particle energies. The n-particle correlated states are labelled by Greek 
letters with the subindex n. For instance as (/32), is a three (two) particle 
correlated state carrying energy W{a-i) {yV{(i'i)). The function A is 

A(m2, J/32; k) = a2^2Y{kj- a2)Y{ki- P2) I ^ 1 , (10) 

y 3 as Oi2 ) 

and 

y(u;«2) = {l + 5{i,3)f'^X{ij-a2). (11) 

The matrix defined in Eq. (|9]) is not hermitian and the dimension may be 
larger than the corresponding shell-model dimension. This is due to the vio- 
lations of the Pauli principle as well as overcounting of states in the CXMSM 
basis. Therefore the direct diagonalization of Eq. dH]) is not convenient. One 
needs to calculate the overlap matrix in order to transform the CXMSM basis 
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into an orthonormal set. In this three-particle case the overlap matrix is 
(0|(c+P+(a2))l3(c;^^(/32)).3|0> = 5.i5a,/3. + $^A(m2,j/32;A;). (12) 

k 

Using the overlap one can transform the matrix determined by Eq. ([9]) 
into a hermitian matrix T which has the right dimension. The diagonalization 
of T provides the three-particle energies. The corresponding wave function 
amplitudes can be readily evaluated to obtain 

|«3) = P+(«3)|0), (13) 
P+(«3) = 5^X(m2;«3)(c+P+(a2))a3, (14) 

iOL2 

where P^{a^) is the three-particle creation operator. 

It has to be pointed out that in cases where the basis is overcomplete 
the amplitudes X are not uniquely defined and, therefore, they do not have 
physical meaning. Instead, the projection of the basis vector upon the corre- 
sponding physical vector, i.e., 

F(za2;«3) = (a3|(c+P+(a2))a3|0>, (15) 

is well defined. 

That the CXMSM wave function amplitudes are usually not well defined 



is no hinder to evaluate the physical quantities. For details see Ref. [2j]. It 
has to be pointed out that this feature does not appear when the basis is 
orthonormal, as in the one-proton three-neutron case to be analyzed below, 
since neither the Pauli principle nor overcounting of states are relevant here. 

The next step in our CXMSM is the evaluation of the three-neutron one- 
proton states. With the basis denoted as 

ba3;a4) = «P+(a3))a4|0), (16) 

where p labels the proton state, 0:3 is as in Eq. f[T^ and are the three- 
neutron one-proton state, the four-particle energies iy(a4) in ^^Li are given 

by 



{W{a,) -Bp- W{as)){a,\{c;P+{as))a,\0) 

5^(pA;;A|V|g/;A)Pi+ ^ {pt; X\V\ql; X)bA 

ql33 {klXa2 ijkl\a2l32 ) 



X (a4|(c+P+(/33))a4|0), 
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where, 



5i = (-lf+^+'=+'X(A:a2;«3)i^(/a2;/33 

X 03/33 P 



p k A 1 J g / A 

Oi2 Ct4 «3 J 1 «2 C(4 



and 



52 = (-ir+^+*+V(jz;a2)r(jA;;/32)X(fca2;«3)F(//32;/33) 

X d,d3/32/33A^ I J' ' Ml.' ^ Mil ' ?V (19) 
\ /32 a4 a3 J \ /32 a4 /33 J \ A; "3 /32 J ^ ^ 

Here p and g label proton states, while i,j,k,l label neutron states. The 
proton-neutron interaction matrix elements {pk; \\V\ql] A) will be discussed 
in the Applications. The wave function amplitudes X and the projected 
quantities F, defined above (Eq. f lTSj) ). have been evaluated in previous 
steps of the CXMSM. Notice that in this case the overlap matrix is the unit 
matrix, i.e., 

(0|(c+P+(a^))t^(c+P+(a3))ajO) = S,,,S^,^, (20) 

The advantage of the MSM in stable nuclei is that one can study the 
infiuence of collective vibrations upon nuclear spectra within the framework 



of the shell model. Thus, in Ref. [2J] the multiple structure of particle- 



vibration coupled states in odd Pb isotopes was analyzed. But the most 
important feature for our purpose is that the CXMSM allows one to choose 
in the basis states a limited number of excitations. This is because in the 
continuum the vast majority of basis states consists of scattering functions. 
These do not affect greatly physically meaningful two-particle states. That 
is, the majority of the two-particle states provided by the CXSM are com- 
plex states which form a part of the continuum background. Only a few of 
those calculated states are relevant, namely the ones that are mainly built 
upon poles. The question of how to evaluate and recognize the physically 
meaningful three-particle states, are addressed in the next Section. 

3. Applications 

In this Section we will apply the CXMSM formalism described above to 
study the nucleus ^^Li. 
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To evaluate the valence shells we will proceed as in Refs. [19|, |22|, |25 
and choose as central field a Woods-Saxon potential with different depths 
for even and odd orbital angular momenta I. The corresponding parameters 
are (in parenthesis for odd Z-values) a= 0.670 fm, rn = 1.27 fm, Vq = 50.0 
(36.9) MeV and Vso=lQ.5 (12.6) MeV. As in Ref. (l9|, we thus found the 
single-particle bound states Osi/2 at -23.280 MeV and 0p3/2 at -2.589 MeV. 
The valence shells are the low lying resonances Opi/2 at (0.195,-0.047) MeV 
and 0(^5/2 at (2.731, -0.545) MeV and the shell 0^/2 at (6.458,-5.003) MeV. 
This cannot be considered a resonance, since it is so wide that rather it is 
a part of the continuum background. Besides, the state lsi/2 appears as an 
antibound state at -0.050 MeV. We thus reproduce the experimental single- 



particle energies as given in Ref. j26|]. We also found other states at higher 



energies, but they do not affect our calculation because they are very high 
and also very wide. We thus include in our Berggren representation only the 
antibound state lsi/2 and the resonances Opi/2 and 0^5/2. 

To define the Berggren single-particle representation we still have to 
choose the integration contour (see Eq. ([T])). 

To include in the representation the antibound lsi/2 state as well as 
the Gamow resonances Opi/2 and 0^5/2 we will use two different contours. 
The number of points on each contour define the energies of the scattering 
functions in the Berggren representation, i.e., the number of basis states cor- 
responding to the continuum background. This number is not uniformity 
distributed, since in segments of the contour which are close to the anti- 
bound state or to a resonance the scattering functions increase strongly. We 
therefore chose the density of points to be larger in those segments. 
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T — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — r 



Figure 1: Contour used to include the antibound state (see, also, Ref. [ill). The points 
Bi denote bound states while A denotes the antiboud state. The points Vi correspond 
to the vertices defining the contour. They have the values yi = (-0.05,0.05) MeV, V2 = (- 
0.1,0) MeV, y3=(0,-0.4) MeV, F4=(0.5,-0.4) MeV, T/5 = (8,-0.4) MeV, F6=(8,0) MeV and 
F7=(10,0) MeV 




Re(E) 



Figure 2: Contour used to include the Gamow resonances represented by the points Gi. 
The vertices are yi=(0,-l) MeV, F2 = (lrl) MeV, y3=(8,-l) MeV, F4=(8,0) MeV and 
V:5 = (10,0) McV. 
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Table 1: Number of Gaussian points in the different segments of the contour of Fig. [T] 
Segment [(0,0)-^!] [V^-V^] [V^-V^] [V^~V,] [V,-V,] [V,-V,] [V^-V,] 
Number 30 30 30 30 30 16 6 



Table 2: Number of Gaussian points in the different segments of the contour of Fig. [2] 
Segment [(0,0)-^] [Vi-V^] [V2-V^] [V^^-V^] [V^-V^] 
Number 30 30 30 8 4 



We include the antibound state by using the contour in Fig. [H The 
number of points in each segment are given in Table [H For the Gamow 
resonances the contour in Fig. [2] is used with the number of Gaussian points 
as in Table |2j 

We have adopted these points after verifying that the results converged 
to their final values. A discussion about the choice of these contours and also 



on the physical meaning of the antibound state can be found in Ref. [20 




Figure 3: Radial function (j){r) corresponding to the single-particle neutron antibound 
state Osi/2 at an energy of -0.050 MeV. 

To explore the extend to which the poles are physically meaningful we 
plotted in Fig. [3] the lsi/2 antibound state. One sees that it extends in 
an increasing rate far out from the nuclear surface, as expected in this halo 
nucleus (the standard value of the radius is here 1.2 x 11^/^=2.7 fm). The 
radial wave function corresponding to the Gamow resonance Opi/2 is shown 
in Fig. m The resonance 0^5/2 has a large and increasing imaginary part at 
relative short distances, as shown in Fig. O 

With the single-particle representation thus defined we proceed to evalu- 
ate the two-particle states. 

3.1. Two-particle states: the nucleus ^^Li. 

We would like to start this Subsection by stressing, once again, that the 
two-particle states that we will investigate here correspond to two-neutron 
excitations in the spectrum of ^^Li, while the corresponding odd proton is 
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2 4 6 8 10 12 14 



r[fm] 



Figure 4: As Fig. [3] for the Gamow resonance Opi/2 at an energy of (0.195,-0.047) MeV. 
The dashed hne is the imaginary part of the wave function. 




Figure 5: As Fig. |3]for the Gamow resonance 0^5/2 at an energy of (2.731,-0.545) MeV. 
The dashed line is the imaginary part of the wave function. 
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an spectator. To avoid confusions in what follows we will call these states 
^^Li(2z/; J), where 2v indicates that the state J, which belongs to the spectum 
of ^^Li, is determined by two-neutron excitations. 

The only state which is measured in ^^Li is its bound ground state, which 
was found to lie at an energy of -0.369 MeV 27|. The corresponding angular 
momentum is 3/2". This spin arises from the odd inert proton, lying deep 
in the spectrum, coupled to two neutrons. The dynamics of the system is 
thus determined by the pairing force acting upon the two neutrons coupled 
to a state 0"*", which behaves as a normal even-even ground state (iol. |22|. 
Besides the energy, this state has been measured to have an angular momen- 
tum contain of about 60% of s-waves and 40% of p- waves, although small 
components of other angular momenta are not excluded [2l|. 

We will perform the calculation of the two-particle states by using the sep- 
arable interaction discussed in Section [2l The strength G\^, corresponding to 
the states with angular momentum A2 and parity (— l)'^^, will be determined 
by fitting the experimental energy of the lowest of these states, as usual. It is 
worthwhile to point out that defines the Hamiltonian and, therefore, is a 
real quantity. The two-particle energies are found by solving the correspond- 
ing dispersion relation while the two-particle wave function components are 
as in Eq. 

With the quantities entering the two-particle TDA equations thus deter- 
mined we evaluated the two-neutron states in ^^Li(2i/; J = 0"*"). We found 
that the angular momentum contain of the ground state wave function is 
46.8% s-states, 49.1% p-states and 4.2% d-states. This is in reasonable agree- 



ment with experiment |21 



The wave function components corresponding to ^^Li(2z/; gs) are strongly 
dependent upon the contour that one uses. However, measurable quantities, 
like the energies and transition probabilities, do not. This is because the 
physical quantities are defined on the real energy axis and, therefore, they 
remain the same when changing contour. But complex states which are part 
of the continuum background do not have any counterpart on the real energy 
axis and the physical quantities for these states acquire different values for 



different contours [20[. We will use this property to determine whether a 
complex state is a meaningful resonance. This is important, since the ground 
state is the only one for which experimental data exists. There might be other 
meaningful states that have not been found yet. This implies that we have 
to evaluate all possible two-particle states which are spanned by our single- 
particle representation. For the purpose of this paper this is an important 
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task, since in the next step of the CXMSM only physically meaningful states 
will be considered as members of the basis. 

To decide whether a calculated state is a meaningful resonance we will 



proceed as in Refs. 20|, |28| and analyze the singlet (S=0) component of the 



two-particle wave function. The corresponding expression for this component 



was given in Eq. (10) of Ref. [28|, but we will show it here again for clarity 
of presentation. For the state a with spin and spin-projection (JM) that 
component is, with standard notation, 

^^ajMim) = [xi/2(l)Xi/2(2)]° 5^ X(a6,«JM)jJfe 

a<b 

X [C{ab,fif2) - {-y^+^''-^C{ba,fif2)] , (21) 

where 

Ciab,m) = Mri)Mr2)i-y'^'^'-'^^'{ f ] [Yuiri)Mr2)]._ 



JM ' 

(22) 

and (f)a{f) is the radial wave function corresponding to the single-particle 
state a (Eq. (^). 

If the two-particle state (aJM) is a meaningful resonance then the wave 
function above should be localized within a region extending not too far 
outside the nuclear surface, and its imaginary part should not be too large 



20[. To study these features it is not necessary to go to all six dimensions 
corresponding to the coordinates ri and 7^2- In fact it is enough to consider 
the coordinate r given by ri = r2 = r = (0,0, r) which corresponds to the 
two particles located at the same point and in the z-direction. For details 
see 28|. We will call this one-dimensional function "^ajAiif)- 

The evaluation of the 0"^ states is a relatively easy task, since in this case 
we have determined the strength Gq+ by fitting the experimental energy of 
^^Li{2i/; gs). With this value of the strength we calculated all the 0"*" states 
and found that the vast majority of them are continuum states which belong 
to the background. The relevant states are those which are built mainly from 
the poles jsf, loj . Although these states may not be in themselves physically 
meaningful resonances, they can infiuence significantly the spectrum of ^^Li. 

We therefore included in the CXMSM basis those two-particle states 
which have at least one pole configuration which, in absolute value, is 0.3 or 
larger. We found that all these states correspond to the poles lsi/2 and Opi/2 
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Figure 6: Radial function 4'(r) corresponding to the two-particle state ^^Li(2z^; 2]^) at an 
energy of (2.300,-0.372)MeV. The dashed line is the imaginary part of the wave function. 

coupled to themselves or to continuum states lying close to threshold. As 
a result, these are narrow states. We found also similar meaningful 1~ and 
2"*" states. Specially important is the state 1~ at (0.084,-0.002) MeV which 
is built by the pole configuration (0si/20pi/2)i- and therefore may be physi- 
cally meaningful. Also the state 2f at (2.300,-0.372) MeV which is practically 
built by the pole configuration (051/20^5/2)2+ can be meaningful. However 
its width (i. e. 0.744 MeV) seems to be too large and therefore we decided to 
analyze the corresponding radial wave function in more detail. For this we 
drew \E'2+(r), as shown in Fig. Ei One sees that the wave function is rather 
localized and that its imaginary part is relatively small as compared to the 
corresponding real part. This is a state which perhaps is at the limit of what 
can be considered a meaningful resonance. Yet, it has an effect on the phys- 
ical three-particle states, as will be seen below. It is worthwhile to point out 
that the width of this state (744 keV) is the escape width. At high energies, 
where the giant resonances lie, most of the width consists of the spreading 
width, i.e., of mixing with particle-hole configurations j29|. However, at the 
low energies of the states that we study this mixing is not relevant. 

An important point for the analysis of the three-particle states to be per- 
formed below, is that the scattering wave functions in the segments [(0, 0) — 
Vi], [Vi — V2] and [V2 — V3] are similar in magnitude to the wave function 
of the antibound state. This is because the segments are very close to the 



antibound state j20|. This is a feature that cannot be avoided, and is due 
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to the attractive character of the pairing force. That is, the lowest single- 
particle configuration in the Berggren basis is V2 ■, with energy — 2e, where 
V2 = (— e, 0). This configuration has to lie above the energy of the two-particle 
correlated state, i.e., it has to be e > u(^^Li{2h'] gs))/2. 

With the states 0"*", 1~ and 2+ thus calculated we proceeded to the cal- 
culation of the three-particle system within the CXMSM. 

3.2. Three-particle states 

Using the single- and two-particle neutron states discussed above, we 
formed all the possible three-particle basis states. Due to the large number of 
scattering states included in the single-particle representation the dimension 
of the three-particle basis is also large. The scattering states are needed in 
order to describe these unstable states. 

With the CXMSM basis thus constructed we evaluated the dynamical 
matrix Eq. ([9]) and the overlap Eq. f|T2|) . With these we formed the sym- 
metric Hamiltonian matrix which we diagonalized to obtain the three-neutron 
states. We found that the lowest state is 1/2+ at (-0.381,+0.023) MeV. That 
is, the energy is real and negative. It is an antibound state, as it is the lsi/2 
state itself. A manifestation of this is that the radial wave function diverges 
at large distances. 

With the three- neutron states thus calculated and the proton state Ops/2 
we constructed the final CXMSM basis to describe the nucleus ^^Li. 

3.3. Four-particle states 

In this case the one-proton three-neutron CXMSM basis is orthonormal 
and the matrix f lT7|) is already the Hamiltonian matrix. All the three-neutron 
quantities appearing in Eq. f[T7|l have been solved in the previous step of the 
CXMSM. Instead, the single-proton energy and the proton-neutron interac- 
tion matrix elements are quantities that we have still to determine. The 
single-proton energy was not considered so far. When we evaluated the 
ground state of ^^Li(2z/) we assumed that its energy was only determined 
by the two-neutron excitations. The value of the energy was obtained by fit- 
ting the strength of the neutron-neutron separable interaction, ignoring any 
effect that the protons may have had. In the analysis of the spectrum of ^^Li, 
including the wave functions, this is irrelevant, since the assumption in those 
calculations was that the protons were only spectators. In other words, the 
effect of the odd proton was only an scaling of the energies. A similar feature 
occurs in the evaluation of the three-neutron one-proton states. As seen in 
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Eq. (fT7|) the value of the proton energy Sp, where p is the proton state 0p3/2, 
only shift the spectrum, but the relative energy between two given states is 
not affected. Again here the effect of the proton degree of freedom regarding 
the single-particle energy is to scale the whole spectrum of ^^Li. 

The determination of Sp is a difficult task due to the energy renormaliza- 
tions that our procedure implies. However, it should be a real number, since 
the proton is a bound state. We will, therefore, not intend to evaluate the ab- 
solute energies and only discuss the ^^Li energies relative to the corresponding 
ground state energy. This is equivalent to take Ep as a real parameter that 
adjust the real part of the ground state energy of ^^Li. The proton-neutron 
interaction matrix elements were taken from empirical effective interactions 



which are determined by fitting experimental data [30|, |31 



With the quantities entering the Hamiltonian matrix f lT7|) determined as 
discussed above, we calculated the energies W{a4) and the wave function 
amplitudes {a4\{CpP~^{a3))a4^\0). Notice that in this case of an orthonormal 
basis the wave function components X and the corresponding projections F 
(Eq. (IT^ ) coincide. 

Due to the presence of the continuum states the dimension of the one- 
proton three-neutron basis is very large. Of all the states calculated within 
this basis we will present only those which are physically meaningful. As we 
have discussed above, the main configurations in such states contain large 
contributions from the poles (relevant configurations). We therefore will 
proceed as in the evaluation of the two-neutron states and consider physically 
meaningful the one-proton three-neutron basis states for which at least the 
amplitude of one relevant configuration is, in absolute value, 0.3 or larger. 

Below we will briefly describe the structure of the states thus calculated 
forming the ^^Li spectrum. 

As seen in Fig. [7] there are five states at low energy predicted by the 
theory. All of them are narrow and mainly built upon configurations that are 
close to the continuum threshold. This is not surprising since the determining 
poles at low energy are the antibound state si/2 and the resonance pi/2, and 
both have real or nearly real energies and are very close to threshold.. 

The ^^Li ground state is 2~ and its energy is real. Since it is neither 
a bound state (its wave function diverges at large distances) neither a res- 
onance (its energy is a real quantity) we conclude that it is an antibound 
state. This is not surprising, since it is mainly built upon the CXMSM basis 
state I [vr(0p3/2)z^(lsi/2 Li(2z/; gs))i/2+]2~), where tt (z/) indicates proton 
(neutron) degree of freedom. The three-neutron component of this state is 
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Figure 7: Experimental level scheme in ^^Li. The three lowest levels are from 
the one at 1.5 MeV is from @. The theoretical results are labelled CXMSM. 



while 



dominated by the state z/(lsi/2) and this induces the state 2]~ to be an- 
tibound. It is worthwhile to point out that the influence of the z/(lsi/2) 
antibound single-particle state is due to its position, lying very close to the 
continuum threshold. Since there is no barrier to trap inside the nucleus the 
neutron moving in this state, the corresponding wave function is very similar 
to the one corresponding to a bound state at the same small and negative en- 



20 . Within the Continuum Shell Model the influence of the antibound 



ergy 

state is taken into account by the scattering states on the real energy axis ly- 



ing close to the continuum threshold [32|. Even in the complex energy plane 



the scattering wave functions lying close to the antibound state are similar 



to each other [20|, l33|]. A discussion on the conveniences and drawbacks of 
using the complex energy plane to describe the antibound state can be found 
in Ref. 
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The question one may ask is how is it possible that in ^^Li three neutrons 
can occuppy the state Is 1/2. The answer to this is that in the continuum this 
state is split in many components. 

The other states at low energy in Fig. [7] are built upon the proton single- 
particle state 0p3/2 coupled to the three-neutron states |lsi/2 Li(2z/; 0"*")), 
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|lsi/2 Li(2z/; 1-)), |0j9i/2 Li(2z/; 0^) and \0pi/2 ®" Li(2z/; 1")). They 
are all narrow states. 

Besides these states there are many levels that lie above 1 MeV which 
are built upon the pole 06/5/2- This may explain the broad resonance seen 
experimentally at about 1.5 MeV. 

The experiment shows two narrow excited states below 0.5 MeV, which 
the calculation predicts to be 1~ and 2+, respectively. 

Most of the calculated levels in Fig. [7] are strongly influenced by the 
continuum states. If we would take only relevant configurations with ampli- 
tudes larger than 0.9 then only two states would appear. If no continuum 
configurations are included then no excited state below 1.1 MeV would be 
found, as can be seen in the shell model calculation of Ref . [7j . 



4. Summary and conclusions 

In this paper we have studied excitations occurring in the continuum part 
of the nuclear spectrum which are at the limit of what can be observed within 
present experimental facilities. These states are very unstable but yet live a 
time long enough to be amenable to be treated within stationary formalisms. 
We have thus adopted the CXSM (shell model in the complex energy plane 
Till) for this purpose. In addition we performed the shell model calculation 



by using the multistep shell model. In this method of solving the shell model 
equations one proceeds in several steps. In each step one constructs building 



blocks to be used in future steps [10|. We applied this formalism to analyze 
the spectrum of ^^Li by assuming that it is determined by one-proton three- 
neutron excitations. First we studied the neutron degrees of freedom by 
profiting from the information that exits on the single-neutron states in ^°Li 
and on two-neutron states in ^^Li. In our case of ^^Li the neutron excitations 
correspond to the motion of three neutrons, partitioned as the one- times two- 
neutron systems. This formalism was applied before, e.g., to study multiplets 
in the lead region 2J| . Finally, the spectrum of ^^Li is calculated by coupling 



the three- neutron system with the 0^3/2 single-proton state. 

We adopted the single-particle energies (i.e., states in ^°Li) as provided by 
experimental data when available or as provided by our calculation. These 
are the antibound state lsi/2 and the resonances Opi/2 and 0^5/2. Besides 
these states (which are poles of the Green function) we also included con- 
tinuum states (Coulomb waves). We found that the physically meaningful 
two-neutron states are mainly built upon poles. With the two-particle states 
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thus obtained we calcuated the three-neutron states and coupled them to the 
0^3/2 protons state to evaluate the spectrum of ^^Li, The ground state energy 
of this nucleus turns out to be real. Since this is neither a bound states (the 
wave function diverge at large distances) nor a resonance (the energy us real) 
we conclude that it is an antibound states. This agrees with a number of 
experiments [s], [g], 0] . 

Besides this antibound state our calculation predicts two low lying narrow 
states which, as seen in Fig. [TJ are also found experimentally. But there are 
also three other calculated levels which have not been observed so far. 

We conclude that the experimentally observed narrow states in ^^Li arise 
as a result of the valence neutrons moving in unstable shells. As in ^^Li, 
these are mainly the antibound state lsi/2 and the narrow resonance Opi/2- 
But in ^^Li continuum states lying close to threshold play also a fundamental 
role. 
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